High energy kink in the single particle spectra of the two-dimensional Hubbard model 
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Employing dynamical cluster quantum Monte Carlo calculations we show that the single particle 
spectral weight A(k,u>) of the one-band two-dimensional Hubbard model displays a high energy 
kink in the quasiparticle dispersion followed by a steep dispersion of a broad peak similar to recent 
ARPES results reported for the cuprates. Based on the agreement between the Monte Carlo results 
and a simple calculation which couples the quasiparticle to spin fluctuations, we conclude that the 
kink and the broad spectral feature in the Hubbard model spectra is due to scattering with damped 
high energy spin fluctuations. 



Introduction. Angle-resolved photoemission spec- 
troscopy (ARPES) has revealed much about the 
cuprates, including the energy scales associated with the 
d-wave gap[l| and a low energy kink presumably asso- 
ciated with strong electron-phonon coupling 2]. Recent 
ARPES experiments have revealed a high-energy (HE) 
kink and a waterfall structure [1,0) 0JH> m which the 
band dispersion broadens and falls abruptly at binding 
energies below « 0.35 eV. The origin of this kink has been 
attributed to a crossover from the quasiparticle (QP) 
to the Mott-Hubbard band0, Q the settlement of spin- 
charge separation Q, or interaction of the q uasi par ticles 
(QP) with spin fluctuation excitations [El [lfl [ll . 12l 13| . 



In this Letter, we study the single particle spectral 
weight A(k, ui) of the one-band 2D Hubbard model with 
near-neighbor hopping t and Coulomb interaction U in 
the regime where U is comparable to the bandwidth 
W = 8t and in the doping range relevant for cuprate 
superconductors. The single-band Hubbard model is be- 
lieved to describe the low-energy physics of the cuprates 
down to energies of 2t below Fermi surface (FS). Sur- 
prisingly, the calculated spectra of the single band model 
are remarkably similar to the experimental ones down to 
binding energies of « At — 5t. They reveal a sharp QP 
feature down to a kink energy Ekink , followed by a steep 
dispersion of a broad waterfall structure. We find that 
these features are accurately captured by a renormalized 
second order (RSO) approximation to the self-energy in 
which the QP couple only to spin fluctuations. A careful 
inspection of the different contributions to the RSO self 
energy shows that the HE kink and the waterfall struc- 
ture is due to the coupling to damped high energy spin 
excitations. 

Formalism. To study the Hubbard Hamiltonian we 
employ the dynamical cluster approximation (DCA)[3, 
[IH, a cluster dynamical mean-field theory which maps 
the original lattice model onto a periodic cluster of size 
N c = L? c embedded in a self-consistent host. Correla- 
tions up to a range L c arc treated explicitly, while those 
at longer length scales are described at the mean-field 
level. With increasing cluster size, the DCA systemati- 



cally interpolates between the single-site dynamical mean 
field (DMFT)[16[ and the exact result. Cluster dynam- 
ical mean field calculations of the Hubbard model are 
found to reproduce many of the features of the cuprates, 
including a Mott gap and strong AF correlations, d- 
wave superconductivity and pseudogap behavior 15j . We 
solve the cluster problem using a quantum Monte Carlo 
(QMC) algorithm 17| and employ the maximum entropy 
method [l8[ to calculate the real frequency spectra. The 
results presented here are obtained from calculations on 
N c = 16 and N c = 24 site clusters for U = 8t. 

Results. The single particle spectral weight of the one- 
band Hubbard model A(k,u>) = —^ImG(k,uj) reveals 
a high energy kink in the QP dispersion. This feature 
is present for a large range of doping values, from the 
undoped system up to w 30% doping and along differ- 
ent cuts of the Fermi surface (FS). In Fig. Q]-a and -b 
we show an intensity map of A(k,oS) along the diagonal 
((0,0) — (it, it)) and center to zone edge ((0,0) — (tt,0)) 
directions at 20% doping. In both cases an intense QP 
peak which cuts the FS can be noticed at small ener- 
gies above Ekink ~ t. At higher binding energies the 
dispersion becomes very steep, the peak broadens and 
decreases in intensity. Ekink is weakly decreasing with 
doping and is weakly dependent on the cut across the FS 
(not shown). These results are in good agreement with 
recent experimental findings 

H II uB, 0, H ■ We find the 
kink position along the diagonal direction to be at a mo- 
mentum larger than (7r/4, 7r/4). The next-nearest and 
next-next-nearest neighbor hoppings t' and t", respec- 
tively, can however modify the position of the HE kink 
in the BZ (not shown). This indicates that (7t/4, 7t/4) 
has no particular relevance for the locus of HE kink in 
the momentum space contrary to some previous sugges- 
tions 0, H, 0| ■ This conclusion is consistent with experi- 
mental results reported for LBCO [j|. 

The HE kink can be inferred from the frequency de- 
pendence of the self-energy S(k, u>). In Fig. [2] we show 
S(k, us) at three different k points along the diagonal di- 
rection. At the kink energy the k dependence of the 
self-energy is weak. Starting from the Fermi energy and 
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FIG. 1: (color) Intensity map of the spectral weight A(k, u>) 
for T = 0.14t, U = 8t and n = 0.80. In (a) and (c), k runs 
along a nodal (0, 0) — (n, n) cut and for (b) and (d) k runs 
along (0,0) — (7T.0). The DCA results shown in (a) and (b) 
were obtained using a 16 site cluster (N c — 16) and the RSO 
results shown in (c) and (d) were obtained using a self-energy 
given by Eq. (1) with U = 0.4(7. A kink followed by a broad 
waterfall feature is noticed below Ekink ~ t. The thin line in 
(a) and (b) indicates the bare dispersion E(k). 
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FIG. 2: (color online) Real part i?e£(k, lu) — fi (thick full) and 
imaginary part — JmE(k,tj) (dashed) of the self-energy at a) 
k = (7r/2,7r/2), b) k = (7r/4,7r/4) and c) k = (0,0). The 
peak in A(k, lu) (thin line) corresponds to the intersection of 
i?eE(k, (J)-n with lu - E(k) (dashed-dotted). The QP is well 
defined down to — lu — Ekink where JmE(k, lu) is small. At 
larger binding energies a waterfall structure develops. 



increasing — u>, i?e£(k, lu) has a negative slope character- 
istic of a QP with an enhanced effective mass. The QP 
is positioned at the intersection of uj — E(k) + (x with 
i?e£(k, lu) and is sharp (see Fig. ®&), a consequence of 
a small ME(k,u). Here -E(k) = — 2t(cosk x + cosk y ) is 
the bare dispersion. This QP feature persists down to an 
energy — lu = Ekink where i?eS(k,w) shows a maximum. 
At larger binding energies i?eS(k, lu) has a positive slope 
which results in the steep dispersion characterizing the 
waterfall region seen in ARPES. The slope increases with 
a finite t' resulting in a steeper waterfall dispersion (not 
shown). However in the waterfall region, /r7iS(k, lu) is 
large yielding a broad and asymmetric feature in A(k, uj) 
(see Fig. [2]-b & c), with the maximum still at the inter- 
section of uj — E (k) + /i with ReYi (k, uj) . The region where 
i?e£(k, lu) has a positive slope spans a large energy range, 
between —t and — 6t, thus characterizing the spectrum 
down to the T point at the bottom of the band (Fig. [J] 
-c). It is interesting to notice that the asymmetry of the 
spectral feature below Ekink in Fig. [2] -b and -c reveals 
the existence of two maxima in A(k,Lu) pushed together. 
In fact these two maxima are much better separated if 
a finite t' is considered, one with a steep dispersion and 
the other with a strongly renormalized one. Similar be- 
havior has been seen in experiment (see Fig. 4 in01). We 
will present results for the Hubbard model with higher 
order hoppings elsewhere. 

The DCA results for the HE kink are different from 
the results of other approaches such as the four-pole 
approximation 1(J U | which considers the scattering of 
the QP in the lower Hubbard band on spin excitations, 
or the DMFT[li|. These studies find a i?e£(k,w) in the 



waterfall region with a slope much larger than one, and 
hence the corresponding A(k,Lu) displays a gap between 
two distinct bands. Here, we find < dReT,(uj)/duj < 1 
near the kink, resulting in a dispersive waterfall feature 
in A(k, lu). 

Since the spin fluctuations are known to be strong 
in the cuprates, a reasonable assumption for the ori- 
gin of the HE kink is the scattering of QP with spin 
excitations (HI. [Tot 11, 12, lj|. In order to understand the 
origin of the HE kink in our results, we therefore test 
how well a simple renormalized second order (RSO) ap- 
proximation to the self-energy given by 
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'(K,iuj) = ^[/ 2 ]T]rG c (K-Q,zc^)x c (Q^) , 

Q " 

(1) 

analytically continued to real frequencies, can describe 
the "exact" DCA self-energy. In Eq. (P), G c (K,iuj) and 
Xc(Q,«i / ) are the interacting cluster DCA Green's func- 
tion and spin susceptibility respectively and U is a renor- 
malized interaction vertex^, [2l|. K and Q are the 
cluster momenta (lij. This approximation neglects the 
frequency and momentum corrections to the interaction 
vertex, and the contributions from the charge and pair- 
ing channels which we find considerably smaller than the 
contribution from the spin channel. By comparing the 
DCA results for A(k, lu) in Fig. [T]-a and b with the corre- 
sponding spectra calculated with the RSO approximation 
shown in Fig. [TJ-c and -d, respectively, one can see that 
the HE kink is well captured by the spin RSO approxi- 
mation. To obtain this agreement we have set U = 0.4(7. 
We also find good agreement between the DCA and RSO 
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results for A(k,ui) at 5% doping with U = 0.3(7 (not 
shown). 

The similarity between the DCA and RSO spectra 
can be deduced from the corresponding self-energies. 
i?e£ °(K, uj) is shown in Fig. [3] -a with dashed lines 
at K = (7r/2,7r/2) and K = (0,0) at 20% doping. Like 
the DCA self-energy (full lines) ReT, RSO (K, uj) shows a 
maximum at lu = —Ekink which causes the kink seen 
in the QP dispersion. The DCA and RSO self-energies 
agree well over the energy range relevant for the HE kink, 
especially in the optimally doped and overdoped regions 
(15% — 30% doping). The agreement is still good at small 
doping as can be seen from Fig. [3]-b, where the 5% dop- 
ing case at K = (tt/A, 7r/4) is shown. However at small 
doping the RSO self-energy gives a smaller Ekink and a 
steeper waterfall dispersion, presumably due to the ne- 
glect of strong coupling effects which become important 
in this region. At positive lu of order of several t the 
RSO self-energy differs from the DCA one, resulting in 
an underestimation of the Mott gap. 

A careful analysis of the different Q contributing to 
the RSO self-energy in Eq. (fT]) shows that the HE kink 
is caused by scattering from high energy spin excita- 
tions. As an example, the red line in Fig. [3] -c is the 
net contribution to the real part of the self-energy at 
K = (tt/2,0) from Xc(Q,v) with Q = (tt,tv/2), Q = 
(0,7r/2), Q = (7r/2,7r/2) and Q = (0,tt)@- It displays 
a maximum at uj = —Ekink, the typical energy associated 
with spin excitations at the magnetic zone boundary (see 
next paragraph), thus yielding the HE kink. The low 
energy spin excitations at the zone center (green) and 
zone corner (blue) do not contribute to the maximum in 
ReT, RSO (K,uj) and therefore are not important for the 
HE kink. The analysis of ReT, RSO (K, lu) at other K val- 
ues leads to the same conclusions. 

Short range spin excitations with a characteristic en- 
ergy 2 J ss 8t 2 /U persist up to relative large doping. The 
magnetic structure factor S(Q,u>) for different values of 
the doping is shown in Fig. 0]-a, -b, -c, -d, & -e at differ- 
ent cluster Q in the BZ. In the undoped system S(Q,lo) 
shows sharp magnon peaks with an energy predicted in 
agreement with the spin-density wave (SDW) approxima- 
tion [23j], as can be seen in Fig. [4]-f where the magnon dis- 
persion along the diagonal direction is shown [24j. With 
increasing doping S(Q,lu) broadens and, in the region 
of the BZ corresponding to high energy spin excitations, 
still retains a well defined maximum at an energy of order 
of w 2 J, as shown in Figs.[4]-a through -d. In this region 
of the BZ the total weight of S(Q,lu) does not change 
much with increasing doping but a significant transfer of 
weight to higher energies can be noticed. For instance the 
peaks in S(Q, uj) at Q = (0, n) and Q = (0, 7r/2) are posi- 
tioned at « 2 J for 0%, 5% and 20% doping. The magnon 
peaks at Q = (7r/2,7r) and Q = (tt/2,tt/2) soften a lit- 
tle with doping, presumably causing the small decrease in 
Ekink with doping. However, we find that the low energy 
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FIG. 3: (color) a) J?eE(k,w) - l uatK= (ir/2,ir/2) (black) 
and K = (0, 0) (red) at 20% doping a) and at K = (tt/4, tt/4) 
and 5% doping b). The full (dashed) lines are the DCA (RSO, 
Eq[T| results, c) Contribution to ReE RSO (K, lo) (black) with 
K = (tt/2,0) from spin excitations with different momentum 
Q. The high energy spin excitations (red) yield a maximum at 
uj = —Ekink- The low energy spin excitations with Q = (0, 0) 
(green) and Q = (n,n) (blue) have a negligible contribution 
to the HE kink. 
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FIG. 4: (color online) a), b), c), d) & e) Doping dependence of 
5(Q, uj) for different Q in the BZ. High energy spin excitations 
persist when the doping is increased, displaying a maximum 
at lo fa 2 J in S(Q,lo). In e) S ((n .n) , uj) is scaled with a factor 
of 0.1. f) Dispersion of the magnon peak along the diagonal 
direction for different dopings. Here U = 8t and T = 0.14t. 



spin excitations around Q = (tt, tt) are more strongly af- 
fected by doping as may be seen in Fig 2] -e. The total 
spectral weight is strongly reduced with doping. At 20% 
doping S((tt, tt), uj) shows a broad peak with a maximum 
at an energy w J. 

Discussions. These results suggest that the HE kink is 
caused by coupling to damped high-energy spin fluctua- 
tions. The dispersive spectral feature in the waterfall re- 



4 



gion is a consequence of i?eS with < dReE(u))/du> < 1. 
This requires scattering on damped excitations with an 
energetically broad spectrum. As seen in Fig(21 the scat- 
tering rate — im£ initially increases with — u). How- 
ever, at higher binding energies JmE passes through a 
maximum and decreases as the phase space for the final 
scattering states decreases. Within the RSO approxima- 
tion this reflects the convolution of the spin-fluctuation 
particle-hole continuum with A(k,oj) and leads, in the 
present case, to a peak in — im£ for lu ks — 3t. The 
structure in the i?eS follows from the Kramers-Kronig 
relation and can be understood a result of energy level 
repulsion. At small values of —to the majority of the 
states in the single-particle-spin-fluctuation convolution 
have energies larger than — u> and give rise to the usual 
QP mass enhancement. However, at larger values of — lo, 
the dominant contribution from these intermediate states 
comes from states with lower energies leading to the de- 
crease in ReYi and driving the dispersion of the spectral 
feature at high binding energies below the bare dispersion 
(thin line in Fig[T]). 

While the main features of HE kink and waterfall can 
be captured with a single-band model with U ~ W, com- 
parison with experiment requires realistic values for the 
Hamiltonian parameters. We already mentioned that 
a next-nearest neighbor hopping t' makes the water- 
fall dispersion steeper, presumably due to sharper spin 
excitations[25]. A t" hopping moves k/ and the locus of 
HE kink on the diagonal direction in BZ towards the T- 
point. We also find that Ekink decreases with increasing 
U presumably due to the reduction of effective J (not 
shown). At high energy, the experimental ARPES in 
cuprates show oxygen valence states in the proximity to 
the r point 0,0], which obviously are not captured with a 
single-band Hubbard model. Moreover other states miss- 
ing in the single band model, such as the d 3z 2_ r 2 states, 
should also be considered when analyzing the experimen- 
tal ARPES spectra below 0.5 eV , as multi-orbital calcu- 
lations for cuprates indicate^. 

The simple renormalized second order ansatz, Eq. [TJ 
seems to provide a good description of the single-particle 
ARPES spectra of the Hubbard model with parameters 
relevant for the cuprates outside the pseudogap regime. 
This suggests that this ansatz could be used to analyze 
experiments where x(<7, w) is measured by neutron scat- 
tering, and used to construct the ARPES spectra. U 
could be fixed by fitting the RSO spectra to the high en- 
ergy kink. Consistency between the measured and con- 
structed spectra would strongly suggest that the HE kink 
in the experimental ARPES spectra may be described 
with a single-band model and is due to the coupling to 
spin fluctuations. As discussed elsewhere 27], a similar 
RSO result also provides an accurate description of the 
pairing interaction of the Hubbard model in the regime 
relevant for the cuprates. Thus the neutron spectra, to- 
gether with the U obtained from the fit above, could be 



used in a simple RSO calculation of the superconducting 
phase diagram. 

Conclusions. By employing DCA calculations we show 
that the single-band Hubbard model captures the HE 
kink structure seen in the cuprates. The kink occurs as 
a crossover from a well defined QP peak to a waterfall 
structure characterized by a broad and asymmetric fea- 
ture with steep dispersion. The structure of the HE kink 
is well captured by a simple renormalized second order 
self-energy which couples the quasiparticle to spin fluc- 
tuations. A careful decomposition of the contributions 
to the RSO self-energy indicate that the HE kink and 
the waterfall structure in the spectrum of the Hubbard 
model is due to the damped high energy spin fluctuation 
continuum. 
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